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SECTION  1 
INTRODUCTION 


The  technique  that  has  evolved  from  this  research  effort  is 
called  HIM  which  stands  for  Hybrid  Iterative  Method.  It  is  the 
third  generation  of  a  technique  that  was  initially  researched 
when  Dr.  G.  A.  Thiele  was  with  the  Ohio  State  University 
ElectroScience  Laboratory.  This  evolution  is  described  below. 

Starting  in  1  979,  Kira  and  Thiele  developed  a  hybrid  AS- MM 
technique  (asymptotic  moment-method  technique)  to  find  the 
induced  currents  on  the  surface  of  scatterers  in  the  intermediate 
frequency  region  [1],  These  induced  currents  consisted  ot  two 
parts:  an  "optics"  current  similar  to  (but  not  the  same 

as  2n  x  ft1);  and  a  "correction"  current  which  was  the  difference 
between  the  optics  current  and  the  true  current.  The  correction 
current  was  obtained  via  the  moment  method. 

In  1983  and  1984,  Kaye,  Murthy  and  Thiele  [2],  [3] 
significantly  improved  on  the  AS-MM  work  by  eliminating  the 
moment  method  regions  thereby  making  the  method  potentially 
applicable  to  large  scatterers  as  well  as  ones  of  modest  size. 

In  the  second  generation  work  by  Kaye,  Murthy  and  Thiele,  the 
magnetic  field  integral  equation  (MFIE)  was  reduced  to  a  sequence 
of  integral  equations  for  both  the  "optics"  currents  and  the 
"correction"  currents,  each  of  which  is  solved  by  iteration.  On 
the  other  hand,  in  Kim  and  Thiele  [1 ] ,  the  "correction"  currents 
were  solved  by  the  moment  method  which  tends  to  limit  the 
electrical  size  of  the  object  that  may  be  considered.  The  second 
generation  work,  then,  was  a  purely  iterative  technique  wherein 
all  the  integral  equations  (tor  both  the  optics  currents  and  the 
correction  currents)  are  of  the  same  form  and  are  similarly 
solved  via  iteration. 
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In  1985  and  1986,  Murthy,  HiLL  and  Thiele  L4J ,  [b]  evolved 

Che  third-generation  technique,  which  is  the  present  technique 
and  which  is  described  in  the  next  section.  This  third- 
generation  technique  is  called  HIM  for  Hybrid  Iterative  Method. 

In  HIM,  we  incorporate  the  correction  currents  into  the  ansatz 
for  the  iterative  scheme  itself  so  that  there  is  only  one  current 
to  iterate.  Thus,  the  iterative  solution  proceeds  by  taking  the 
previously  iterated  current  and  using  it  as  the  "new"  current  to 
be  improved  upon.  Typically,  ten  iterations  is  sufficient  to 
obtain  a  converged  correct  current.  The  second-  and  third- 
generation  techniques  were  mostly  supported  by  RADC/Hanscom . 

The  salient  features  of  the  hybrid-iterative  method  (HIM) 
may  then  be  stated  in  the  following  way.  HIM  is  an  iterative 
algorithm  to  solve  MFIE  (magnetic  field  integral  equation)  tor 
induced  surface  currents  on  a  scatterer.  Every  iterative  scheme 
must  be  started  with  an  initial  guess.  The  better  the  initial 
guess,  the  faster  the  method  converges.  This  initial  guess  must 
be  obtained  from  a  knowledge  of  the  physics  of  the  scattering 
process  to  the  extent  it  is  known.  Thus,  in  hIM,  the  physics  ot 
the  scattering  process  is  incorporated  into  the  solution  process 
at  the  outset  to  the  extent  it  is  possible,  in  the  form  ot  an 
initial  guess  for  the  current  on  the  shadow  side.  This  results, 
directly,  in  the  reduction  of  CPU  time.  It  is  to  be  noted  that 
the  iteration  method  [2],  [3]  is  a  special  case  of  HIM  in  that  it 
takes  the  shadow-side  current  to  be  zero.  For  some  scattering 
geometries,  such  a  crude  initial  guess  may  not  be  good  enough  to 
warrant  a  converged  solution.  In  HIM,  the  solution  always 
converges  when  the  initial  guess  is  other  than  zero. 

In  Section  II,  we  develop  the  integral  equations  tor  HIM  and 
discuss  their  solution  by  iteration.  In  Section  III,  we  present 
theory  relevant  to  the  initial  estimates  ot  the  shadow-side 
currents  tor  both  smooth  and  edged  bodies.  Section  IV  presents 
calculated  results  using  the  theory  developed  in  section  II. 

Both  edged  bodies,  exemplified  by  a  square  cylinder  and  smooth 
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bodies,  exemplified  by  circular  and  elliptic  cylinders  are 
considered.  Section  IV  also  presents  our  results  obtained  on  tne 
Cray  X-MP.  In  addition,  ogival  cylinders  which  involve  Doth 
edge-  and  smooth-surface  diffraction,  and  hence  constitute  a 
stricter  test  of  the  method,  are  also  considered.  Section  V 
contains  the  summary  of  the  technique  and  points  out  further 
avenues  of  research. 
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SECTION  2 

GENERAL  HIM  THEORY 

2.1  MAGNETIC  FIELD  INTEGRAL  EQUATION 

Consider  a  perfectly  conducting  body  illuminated  by  a  plane 
wave.  The  induced  surface  current  density  on  the  body  at  an 
observation  point  P,  as  shown  in  Figure  1 ,  can  be  computed  from 
the  MFIE: 

J(ft)  =  2n  x  ft  (ft)  +  2n  x  L[J]  (1) 

where  the  operator  L  is  defined  by 

L[J]  =  /  3(ft')  x  $'G(r)  ds'  .  (2) 

Z 

The  surface  of  the  body  is  denoted  by  Z  and  ft  and  ft'  are  the 
observation  and  source  point  position  vectors,  respectively, 
on  Z  and  r  =  |ft  -  ft' |  .  The  prime  on  the  gradient  operator 

indicates  that  the  differentiation  is  performed  on  the  source 
coordinates.  The  bar  through  the  integral  sign  is  used  co  denote 
the  principle  value  integral  over  Z  .  ftL  is  the  incident  magnetic 
field  vector  and  n  is  the  outward  unit  normal  to  the  surface 
at  ft.  G(r)  is  the  tree-space  Green's  function  given  by 

e-jBr 

G(r)  =  - -  (3) 

4  it  r 

for  the  three-dimensional  problem  and  by 

G(r)  =  J-  H<2> (8r)  (4) 

4j 

for  the  two-dimensional  problem.  H^^(Sr)  is  the  zero-order 
Hankel  function  of  the  second  kind  and  S  is  the  free  space 
propagation  constant.  The  time  dependence  is  taken  to 
be  exp(jcut)  and  is  suppressed  throughout. 
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2.2  CURRENTS 


The  first  step  in  solving  the  MFIt-  tor  the  induced  surtace 
current  density  J(R)  is  to  divide  the  total  surface  £  of  tne  oooy 
into  £*  and  Is,  which  represents  the  lit  and  shadowed  regions, 
respectively.  The  dividing  line  between  these  two  regions  is  tne 
geometrical  optics  shadow  boundary  which  is  defined  as  the  locus 
of  points  satisfying  u  n  =  0,  where  u  is  the  unit  vector  in  the 
direction  of  propagation  of  the  incident  field.  Then  3(R)  can  be 
expressed  as  follows: 

J(ft)  =  63i(R)  +  (1  -  <5):fS(&)  (5) 

where  6=1  if  &  e  E*  and  6=0  if  fi  £  Es.  Eq.  (1)  can  now  be 
rewritten  as 

3(R)  =  2n  x  fL(R)  +  2n  x  L[  ]  +  2n  x  L[JS]  .  (6) 

In  the  notation  used  in  Eq .  (6)  and  in  subsequent  equations,  the 
surface  of  integration  is  identified  by  the  superscript  on  the 
current  density.  We  now  write  out  Eq .  (6)  explicitly 
for  H  c  and  Re  E s : 


=  2n 

x 

+  2n  x 

ElJ*] 

+  Zn  x  L[ 3s  ] 

(7) 

3S(R) 

=  2n 

x  a.  (ft) 

+  2n  x 

L[Jl] 

+  2n  x  L[jsj 

(«) 

The  problem  of  solving  for  the  induced  surface  current 
density  has  now  been  transformed  into  one  of  solving  Eqs.  (7)  and 
(8)  for  3*  and  3s.  We  compute  these  currents  in  sequential 
fashion  in  the  following  way. 

First,  let  Jg  be  an  estimate  of  shadow  current.  For  edged 
bodies  like  the  square  cylinder,  this  current  is  obtained  from 
the  wedge  diffraction  theory.  For  smooth  bodies  like  circular  or 
elliptic  cylinders,  this  current  is  obtained  using  Fock  theory. 
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A  detailed  discussion  of  determining  will  be  considered  in 
Section  III.  Thus, 


(9) 


Substituting  tq .  (9)  in  tq.  (7)  and  noting  that  tne  current  on 
the  lit  side  would  now  be  an  approximation  to  the  true  current 
and  denoting  this  zeroth  order  approximate  current  by  3*, 

}*(ft)  =  2n  x  ftL(&)  +  2n  x  L[Jjj  +  2n  x  J  (10) 

This  current  is  a  significant  improvement  over  the  classical 
physical  optics  (PO)  current  since  it  takes  into  account  mutual 
interaction  of  current  on  both  lit-  and  shadow-regions  in 
addition  to  the  geometrical  optics  (GO)  field. 

A  shadow-side  zeroth  order  approximate  current,  J®,  is  now 
defined  by  substituting  J*  =  Jj*  in  (8): 

J*(ft)  =  2n  x  ft.  (ft)  +  2n  x  L[Jj]  +  2n  x  L[(J®]  (11) 

In  this  expression,  the  main  contribution  of  the  second  term  on 
the  right-hand  side  will  be  approximately  -2n  x  ft^(ft)  to  cancel 
the  first  term.  The  3®  obtained  from  tq.  (11)  would  be  closer  to 
the  true  current  than  our  initial  guess  J ® .  Indeed,  the  only 
inaccuracy  incurred  being  that  through  the  approximate  value 
taken  for  the  lit-side  current. 

Now  that  an  approximation  to  the  shadow-side  current  better 
s 

than  Jg.  is  available,  the  next  higher-order  current  on  the  lit 
side  may  be  obtained.  This  lit-side  current  may  be  used  to 
further  improve  the  shadow-side  current.  Thus,  higher-order 
currents  may  be  obtained  until  the  desired  accuracy  is  achieved. 
The  integral  equations  for  the  nth  order  currents  may  be  defined 
as  follows: 
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J*(ft)  =  2n  x  +  L[j*]  +  ] 

3®  (ft)  =  2n  x  ft.  (ft)  +  L[JJ]  +  L[J®] 


(12) 


(13) 


n  =  1  ,  2  .  .  . 

is  the  solution  of  Eq.  (11). 

In  the  iterative  technique  described  in  an  earlier  paper 
[2],  the  "optics"  currents  were  improved  upon  by  the  addition  ot 
correction  currents.  The  technique  described  here  is 
mathematically  equivalent  to  the  correction  current  ansatz  or 
that  iterative  technique  and  at  the  same  time  is  simpler  from  the 
programming  point  of  view. 

2.3  CONVERGENCE  CRITERION 

We  now  give  a  criterion  for  the  convergence  of  the  sequence 
of  integral  equations.  This  criterion  is  based  upon  how 
accurately  our  solution  satisfies  the  MFIE.  Using  Equation  (5), 
we  define, 

Jn(ft)  =  «  (ft)  +0  -  «)  3nS(ft) 

Then  the  current,  Jn(ft)  ,  is  obtained  Dy 
MFIE.  That  is, 

Jn(ft)  =  2n  x  Hi(S)  +  2n  x  L[Jn(K)]  (15) 

If  Jn(ft)^were  to  be  exact,  then  Jn(£)  =  Jn(K).  The  "ditference" 

between  J  (ft)  and  .1  (ft)  is  a  measure  of  the  degree  of  accuracy  or 

the  solution  after  n-orders  of  iteration.  Noting  that  the 

currents  are  specified  as  a  set  of  N  points,  the  average 

★ 

error,  e,  between  Jn  and  Jn  is  defined  to  be, 


(14) 

substituting  J  (ft)  in 
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]_ 

N 


N 

E 

i=1 


1/2 


(16) 


€ 

n 


If  the  error  e  is  less  than  a  preselected  value,  the  sequence 
may  be  considered  to  have  converged. 

Thus,  the  hybrid-iterative  method  (HIM)  computes  the 
currents  on  the  lit  and  shadow  regions  in  sequential  fashion 
starting  with  an  initial  estimate  tor  the  shadow  current.  All 
the  integral  equations  are  of  the  same  form  and  are,  in  tact, 
Fredholm  integral  equations  of  the  second  kind.  Hence,  all  ttiese 
integrals  may  be  solved  by  iteration  as  discussed  in  [2J,  13]. 
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SECTION  3 
INITIAL  CURRENTS 


An  initial  estimate  of  the  current  on  a  scatterer  with 
surface  discontinuities  like  a  square  cylinder  is  obtained  from 
the  knowledge  of  currents  on  a  wedge.  For  a  scatterer  with 
smooth  surface,  Fock  theory  may  be  used  to  obtain  the  initial 
estimate  of  the  shadow-side  current. 


3.1  SHADOW-SIDE  CURRENT  ON  A  WEDGE 


Consider  the  wedge  shown  in  Figure  2  illuminated  by  a  TE- 
plane  wave.  When  the  angle  of  incidence  is  such  that  the 
shadow  boundary  is  not  "close"  to  tace  B ,  the  current  on  race  a 
is  given  by  [ 5  J 


s 

E 


where 


H 

4  —  K  (x  ) 
N  ’  B 


o  <  p  <  <*> 

o  <;  4>i  (n-1  ) w 


— 

10 


(17) 


.  2 


■j  c2 


K_(x)  =  /  1  eJX  /  e"Jk-  dt 

IT  X 

. 

Xg  =  /  2SPg  (cos  —  +  cos  — )/sin  — 

n  n  n 

H0  =  amplitude  of  the  incident  magnetic  field 


p  =  distance  from  the  edge 
n  =  (2  -  — ) 

IT 

a  =  interior  wedge  angle. 


When  the  angle  of  incidence  is  such  that  a  shadow  boundary  is 
close  to  Face  B,  the  current  is  given  by, 
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(18) 


Thus,  Equations  (17)  and  (18)  express  the  currents  on  the 
shadow  side  of  a  wedge  in  closed  form.  These  expressions  involve 
the  well-known  modified  Fresnel  tunctions,  K._(x)  [bj,  and  are 
easy  to  compute.  Figure  3  presents  shadow-side  currents  on  wedge 
computed  using  the  above  expressions  as  well  as  the  exact 
solution  obtained  from  the  eigenfunction  solution.  Note  that  the 
agreement  is  excellent. 

These  currents  on  the  wedge  may  readily  be  used  to  obtain  an 
initial  estimate  of  the  shadow-side  current  on  a  square 
cylinder.  Consider  the  geometry  of  a  square  cylinder  shown  in 
Figure  4.  Faces  2  and  3  are  in  the  shadow  region.  Faces  1  and  2 
constitute  a  wedge  illuminated  with  the  angle  of  incidence 
being  g.  The  current  on  Face  2  may  be  approximated  by  the 
current  on  a  corresponding  wedge  with  a  =  ir/2  and  In 

1  1  ,  D 

a  similar  fashion,  the  current  on  Face  3  also  may  be 
approximated . 


3.2  SHADOW- SIDE  CURRENT  ON  A  SMOOTH  bODY 

An  estimate  for  current  on  the  shadow  side  of  a  smooth 
convex  cylinder  may  be  obtained  by  using  Fock.  theory  [4  j  . 
Referring  to  Figure  5,  Oi  ancl  0'2  are  Roi-hts  of  grazing  incidence 
and  Oq  is  the  observation  point  where  current  is  to  be  found. 

The  field  at  Qy  is  due  to  two  creeping  waves  launched  from 
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Exact  Solution 
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_  fcxact  Solution 

x  x  x  x  Closed  to rm  Solution 


Q2  and  traveling  the  distance  of  and  t2  on  the  surtace, 
respectively.  Then,  for  TE  polarization,  using  the  generalized 
Fock  theory  [7],  [8],  one  obtains 

*  *  -jet,  .  -jet2 

(Q0)  =  ^(Q^e  g(C  1  )  +  fii(Q2)e  Z  g(c2)  (19) 


where 


Qn 

2-/  ”<C> 

’  Ql  .2  PgCf) 

(20) 

m(t’)  =  [SPg(t')/2]1/3 

(21) 

% 

t  2  =  /  df 

1  •  ^  M 

(22) 

g(C) 


I  e-j^c  -dc 
c  w^(t) 


(28) 


P 


w. 


(t')  is  the  radius  of  curvature  at  t'  and  the  Airy  function 
S 

(t)  and  the  contour  "c"  are  defined  in  [4]. 


Once  the  magnetic  field  is  found,  the  current  on  the  surface 
is  determined  from 


3  =  n  x  ft  . 


(24) 


For  the  specific  case  of  a  circular  cylinder  excited  by  a  TE- 
plane  wave,  the  initial  estimate  for  the  shadow-side  current,  3®, 
is  given  By 


3^  =  n  x  ft 
E 


(25) 


1  7 


with  ft  being  given  by 


UbJ 


+  -  -jet.  -jec2 

H  =  zHQ{e  g  ( C  ,, )  +  e  g(?z)) 


j- £a_^  1  /3  ^1_ 
^  2  a 

( ea.)  1  /3  4 
2  a 


a  =  radius  of  circular  cylinder. 


(27) 


Using  Equations  (25)  to  (27),  Fock  current  in  the  shadow 
side  has  been  computed  for  a  circular  cylinder  of  radius  3.2i  and 
is  shown  in  Figure  6.  Notice  that  expressions  given  here  are 
valid  only  for  TE  polarization.  Similar  expressions  can  De 
easily  obtained  for  TM  polarization. 
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Method  ot  Moments 


SECTION  4 

NUMERICAL  EXAMPLES 


In  this  section,  we  illustrate  the  technique  described  in 
earlier  sections.  We  deal  with  scatterers  with  surface 
discontinuities  as  well  as  smooth  surfaces.  Thus,  1)  scatterers 
with  edges  are  exemplified  by  a  square  cylinder;  2)  scatterers 
with  smooth  surfaces  are  exemplified  by  circular  and  elliptical 
cylinders;  and  3)  scatterers  with  a  smooth  surface  joined  to  an 
edge  are  exemplified  by  the  ogival  cylinder. 

Before  the  theory  presented  in  earlier  sections  can  be  used 
to  compute  the  induced  currents,  some  numerical  aspects  must  be 
considered.  For  scatterers  with  edges,  special  attention  must  be 
paid  to  evaluating  the  integrals  for  current  very  close  to  the 
edge.  This  analysis  has  been  dealt  with  elsewhere  [2],  [3]  and 

will  not  be  repeated  here.  For  scatterers  with  smooth  surfaces, 
it  is  necessary  to  evaluate  integrals  wherein  the  integration  and 
observation  points  may  coincide.  In  such  caser ,  the  Green's 
function  becomes  singular  and  one  must  ousider  the  principle 
value  ot  the  integral.  For  two-dimensional  geometries,  however, 
it  can  be  shown  [2],  13]  that  the  integrands  nave  a  finite  limit. 

The  following  results  will  be  tor  zeroth  order,  tirst-order 
and  second-order  currents.  It  is  necessary  to  understand  what 
this  means.  The  zeroth-order  currents  are  defined  by  tqs.  (10) 
and  (11).  The  zeroth-order  currents  are  themselves  obtained  by 
iteration,  say  ten  iterations.  Once  a  converged  zeroth-order 
current  has  been  obtained  according  to  the  criterion  of  Eq .  (16), 
the  zeroth-order  current,  is  used  in  Eqs.  (12)  and  (13)  to  obtain 
the  first-order  currents  via  iteration,  etc. 

4.1  SQUARE  CYLINDER 

Two  examples  are  presented  here.  In  both  cases,  the  size  ot 
the  cylinder  is  the  same,  viz,  w  =  3.7X.  However,  the  angles  ot 
incidence,  are  different.  The  case,  $  =  115°,  has  been 

2U 


dealt  with  in  an  earlier  report  [3].  We  study  the  same  case  tiere 
using  our  hybrid-iterative  method  (HIM).  It  was  necessary  in  [3] 
to  compute  tour  orders  ot  correction  currents  to  obtain  a 
sufficiently  accurate  result.  With  the  HIM,  zeroth  order 
currents  themselves  give  a  sufficiently  accurate  solution,  we 
also  present  the  first-order  currents  for  comparison.  Figures  7 
to  11  show  the  results  for  w  =  3.7X  and  <p  ^  =  115°. 

We  also  consider  the  case  when  ^  =  95°.  Note  that  in  this 
case,  Face  2  of  the  square  cylinder  is  almost,  but  not  quite,  in 
the  visible  region.  As  can  be  seen  from  the  final  result,  the 
current  on  this  face  is  quite  appreciable  and  may  not  be 
neglected  initially  as  was  done  in  the  iterative  technique  [2], 
[3].  Figures  12  to  16  show  results  corresponding  to  this  case. 

In  this  more  demanding  case  it  was  necessary  to  compute  up  to  the 
second-order  currents. 

4.2  SCATTERERS  WITH  SMOOTH  SURFACE 

In  [2],  Kaye,  Murthy  and  Thiele  used  the  iterative  method 
(IM)  to  compute  the  surface  currents  on  a  circular  cylinder  ot 
radius  3.2X.  These  computations  include  the  first-order 
correction  current  and  are  characterized  oy  a  ripple.  For  the 
same  cylinder,  we  computed  induced  currents  using  the  HIM.  These 
results  are  presented  in  Figures  17  to  19.  Note  that  there  is  no 
ripple  and  that  the  agreement  between  our  results  and  the  exact 
eigenfn.net  ion  solution  is  excellent.  Furthermore,  we  needed  to 
compute  only  up  to  the  zeroth-order  currents.  After  zeroth 
iteration,  HIM  yielded  an  error  of  eQ  (HIM)  =  0.0018.  The  VAX 
CPU  time  taken  is  2  minutes  and  18  seconds.  For  the  same  case, 
with  iteration  method,  eQ  (IM)  =  0.054  and  e1  (IM)  =  0.047. 

Thus,  even  after  the  first-order  calculation,  the  average  error 
for  IM  is  significantly  larger  than  that  for  HIM.  The  CPU  time 
required  for  IM  up  to  first-order  calculations  is  three  minutes 
and  45  seconds.  All  computations  were  carried  out  on  a  VAX 
11/780. 
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Figure  19.  Bistatic  echo  width  ot  a  circular  cylinder 
Order  -  0,  Radius  3 .2  A,  *.  =  90° 


In  order  Co  demonstrate  that  our  technique  works  well  tor 
"small"  bodies,  we  also  considered  a  0.21  radius  circular 
cylinder  in  [9].  We  computed  both  zeroth-order  and  tirst-order 
currents  with  excellent  agreement  between  our  results  and  the 
exact  results  obtained  using  the  eigenlunct ion  solution. 

The  computer  code  we  have  written  is  general  and  may  be  used 
to  compute  the  induced  currents  on  scatterers  with  an  arbitrary 
convex  cross-section.  To  demonstrate  this  generality,  we  have 
also  considered  an  elliptic  cylinder  as  shown  in  Figure  2(J  with 
serai-major  and  serai-minor  axes  being  1.51  and  1.01, 
respectively.  Our  results  are  compared  with  those  of  the  method 
of  moments  (MM) .  A  result  is  shown  in  Figure  21 .  Note  the 
excellent  agreement  between  the  results  obtained  with  the  MM  and 
and  with  the  hybrid-iterative  technique. 

4.3  SCATTERERS  WITH  A  SMOOTH  SURFACE  JOINED  TO  AN  EDGE 

The  ogival  cylinder  is  an  example  of  a  smooth  surface  joined 
to  an  edge.  To  illustrate  the  application  of  HIM  to  the  ogival 
cylinder, we  will  consider  a  3A  long  (edge-to-edge)  ogival 
cylinder  with  an  80°  interior  "wedge"  angle  at  each  edge.  Figure 
22  shows  the  current  on  one  side  of  the  ogival  cylinder  tor  a 
plane  wave  incident  at  9U° .  The  current  is  of  zeroth  order  and 
the  agreement  is  seen  to  be  excellent.  The  traveling  wave,  tor 
which  the  ogive  is  well-known,  would  be  predicted  by  the  current 
in  Figure  22.  Figure  23  shows  the  ogival  cylinder  illuminated 
broadside  (i.e.,  at  U° ) .  This  is  a  more  difficult  case  than 
with  90°  (edge  on)  incidence.  Thus,  second-order  currents  are 
required  to  obtain  excellent  agreement. 

4.4  VERY  LARGE  SCATTERERS 

For  the  purpose  of  adapting  the  HIM  to  creat  scatterers  chat 
are  electrically  large,  the  surface  current  is  considered  to 
consist  of  two  components:  the  physical  optics  current  and  the 
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Figure  20. 


Geometry  for  an  elliptic  cylinder  with  semi- 
tnajor  axis  "a"  ana  semi-minor  "b”.  .  «  yyo 
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igivai  cylinder  results  tor  a  JA  long  au°  ogive  wild 
dge-on  incidence.  Current  is  zeroth  order  only, 
xceilent  agreement  between  HIM  and  riM. 
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nonuniforra  current.  The  physical  optics  current  is  readily 
determined.  Since  the  nonuniform  current  is  induced  by  the 
diffracted  field  at  che  surface,  it  decays  relatively  quickly 
from  the  point  of  diffraction.  Thus,  though  the  scatterer  is 
large,  the  region  over  which  the  nonuniform  current  is  unknown  is 
relatively  small.  This  is  the  basic  principle  behind  the 
reformulation  of  HIM  for  large  bodies.  Additionally,  the 
integrals  involved  are  in  a  form  suitable  tor  evaluation  oy  the 
stationary-phase  method;  this  form  is  used  wherever  feasible. 
Thus,  numerical  integration  is  used  to  the  least  extent 
necessary.  We  give  the  specific  formulation  tor  a  circular 
cylinder  in  [10],  though  the  method  itself  is  more  general  than 
presented  there. 

Figures  24  and  25  illustrate  the  foregoing  principles  by 
presenting  a  comparison  of  the  current  between  the  exact 
eigenfunction  solutions  and  the  calculated  results  obtained  by 
use  of  this  reformulated  scheme  for  the  circular  cylinder  with 
radii  equal  to  25  and  80  wavelengths.  The  stationary-phase 
evaluation,  in  the  result  shown,  only  includes  the  endpoint 
contribution  because  the  stationary  points  are  outside  the 
integration  path.  The  CPU  time  required  to  run  the  result  of 
this  25-wavelength-radius  cylinder  on  a  VAX  11/78U  computer  is 
approximately  eight  minutes;  while  the  CPU  time  required  tor  a 
cylinder  with  an  80-wavelength  radius  is  approximately  40 
minutes . 

The  large  body  formulation  tor  tilM  is  adaptable  tor  use  on  a 
super  computer  like  th  Cray  X-MP.  On  May  19,  20,  and  21,  1986, 
we  ran  our  program  for  the  large  cylinder  on  the  Cray  X-ivip  at  ttie 
Cray  facility  in  Minneapolis  under  the  sponsorship  of  Cray 
Research.  Some  rewrite  of  our  program  was  necessary  to  permit 
che  compiler  to  vectorize  the  code.  The  results  are  presented  in 
the  table  which  follows.  Note  that  for  the  80X  cylinder  the  VAX 
time  was  about  2,400  seconds  and  that  on  the  Cray  X-MP,  after 
rewrite  and  vectori zat ion ,  the  CPU  time  was  but  6.17  seconds. 
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RCS 

COMPUTER  TIMES  FOR  LARGE  CIRCULAR  CYLINDERS 


RADIUS  (X) 

~C(  X  ) 

VAX* 

CRAY  X-MP* 

KATIU 

10 

63 

60 

0.37 

162 

25 

157 

480 

80 

500 

2.400 

6.17 

3ay 

1  20 

750 

10.5 

160 

1 ,000 

15.3 

200 

1  ,250 

20.5 

240 

1 ,500 

26.0 

300 

1  ,880 

34.8 

400 

2,500 

50.6 

440 

2,760 

56-.  9 

TIME  IN 

SECONDS 
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While  this  ratio  of  389  to  1  was  not  the  largest  ratio  the  Cray 
Applications  Group  had  ever  seen,  it  was  in  the  class  of  very 
large  ratios  that  one  can  achieve  in  going  from  a  VAX  to  a 
Cray.  They  viewed  a  ratio  of  100  or  over  as  very  good. 

Note  Chat  the  small  CPU  times  indicate  that  the  hIM  could  be 
used  as  an  interactive  design  tool.  For  example,  one  could 
design  any  wing  profile  and  any  radar  frequency  interactively 
using  this  HIM  approach. 

Note  that  we  used  here  the  (classical)  circular  cylinder 
shape  because  we  have  an  exact  solution  against  wnich  to 
compare.  Cylinders  of  other  cross  section  shapes  would  have 
similar  running  times. 

4.5  PERFECTLY  CONDUCTING  CUBE 

Figure  26  shows  the  geometry  for  a  perfectly  conducting 
cube.  Figures  27-32  show  the  H-plane  and  E-plane  bistatic 
scattering  pattern  for  cubes  that  are  0.75X,  1 . 5 X  and  3. OX  on  a 
side,  respectively.  In  all  cases,  agreement  with  experimental 
data  furnished  by  RADC/Hanscom  is  seen  to  be  quite  good.  Also 
included  are  figures  33-36,  which  describes  the  surface  current 
densities  on  the  cube  faces.  Figure  33  describes  the  magnitude 
of  the  z-component  of  the  initial  current  on  the  side  face,  while 
figure  34  shows  that  same  current  aften  imposition  of  the 
M.F.I.E.  It  is  seen  that  the  composite  effects  of  slope 
diffraction  and  vertex  diffraction,  previously  ignored,  have  been 
introduced.  Figures  35  and  36  describe  current  components 
initially  assumed  to  be  zero.  After  imposing  the  M.F.I.E.,  it  is 
seen  that  a  much  better,  non-zero  approximation  results. 

4.6  TM  CASE  FOR  SQUARE  CYLINDER 

Our  results  in  Section  4.1  tor  the  square  cylinder  were  tor 
the  TE  case.  Clearly,  to  do  the  cube  in  Section  4.5  it  was 
necessary  to  be  able  to  calculate  TM  currents  as  well  as  TE 
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currents.  Here,  in  this  section,  we  will  briefly  present  some 
results  for  TM  currents  on  a  square  cylinder. 

Figure  37  shows  the  magnitude  of  the  fourth-order  current  on 
a  square  cylinder  3.193X  on  a  side.  The  incidence  angle  is 

=  135°  (see  Figure  4).  Figure  38  shows  the  phase.  Figures  39 
and  40  show  similar  results  except  that  the  incidence  angle  is 
120°. 
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Figure  i 7.  Magnitude  ut  current  on  a  square  cylinder.  Order 

/* ,  w  -  :j.  i<m  .  + .  *  m° . 
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Figure  33.  Phase  ut  current  on  a  .square  cylinder.  Order 
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D I STANCE ( WAVELENGTHS ) 

Figure  39  .  Magnitude  ot  current  on  a  square  cyLinder.  Order 
9,  w  =  j  .  1  9  3  X  d>  .  =  1  ZUu  . 


D I STANCE ( WAVELENGTHS ) 

Figure  40.  Fhase  of  current  on  a  square  cyiinder. 
w  *  3. 193X  .  • .  *  120° . 
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SECTION  5 


SUMMARY 

A  new  iterative  method  has  been  presented  tor  computing  the 
current  induced  by  pLane  wave  excitation  on  conducting  bodies  of 
arbitrary  shape.  In  this  method, the  scattering  body  is  divided 
into  lit-  and  shadow-side  regions  separated  by  the  geometric 
optics  boundary.  An  MFIE  is  written  for  each  region.  Each  MFIE 
is  solved  by  iteration  (i.e.,  method  of  successive 
approximations).  In  order  to  accomplish  this,  it  is  often 
necessary  to  have  an  initial  estimate  of  the  shadow-side 
current.  This  estimate  is  obtained  from  Known  physical 
scattering  phenomena  (e.g.,  edge  ditrraction  theory  or  Fock 
theory).  hence,  the  method  is  called  a  hybrid-iterative  method 
(HIM)  . 

The  KIM  is  well-suited  tor  use  on  a  super  computer  like  the 
Cray  X-MP.  Results  were  shown  for  2-dimensional  objects  whose 
size  varied  from  0.6X  to  44UA ,  all  with  CPU  times  under  one 
minute,  on  the  Cray.  With  such  speed,  it  is  possible  to  have  a 
computer  code  with  which  interactive  Aero/RCS  tradeoffs  can  be 
evaluated . 

The  HIM  method  is  a  current-based  approach  to  electro¬ 
magnetic  scattering.  In  a  current-based  method,  the  scattered 
field  is  obtained  from  Che  induced  currents.  This  can  be 
contrasted  with  the  field-based  approach  (e.g.,  GTD)  wherein  the 
scattered  field  is  obtained  from  the  incident  field  via  ray 
tracing  without  Che  necessity  of  explicitly  determining  the 
current.  For  modeling  large  comp  lex  geometries,  it  is  believed 
chat  the  current-based  approach  is  better  chan  the  field  based* 
approach  because  one  does  not  need  to  trace,  perhaps,  hundreds  on 
rays  and  one  does  not  need  to  describe  geometrical  surfaces  to 
high  precision. 

The  HIM  approach  to  electromagnetic  scattering  has  evolved 
from  eight  years  of  research  sponsored  primarily  by  KADC/Hanscom 
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with  early  support  trora  ONK  as  weLl  as  lesser  support  trom  wpaI'b 
and  Rockwell  Sciences  Center.  The  merits  of  HIM  have  been  well- 
demonstrated  for  perfectly  conducting  uncoated  two-dimensional 
scatterers.  A  number  of  fundamental  problems  remain  to  be 
researched,  however,  before  the  technique  can  be  used  by  the  KCS 
community  to  model  coated/partially  coated,  three-dimensional 
shapes . 
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l 

1  MISSION  2 

$  of  l 

§  Rome  Air  Development  Center  j 

*  RADC  plans  and  executes  research,  development,  test  and  selected  ^ 

acquisition  programs  in  support  of  Command,  Control,  Communications  3> 

&  and  Intelligence  (C3!)  activities.  Technical  and  engineering  support  within  X 

*•  areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other  ( 

3  ESD  elements  to  perform  effective  acquisition  of  C3 1  systems.  The  areas  „ 

g.  of  technical  competence  include  communications,  command  and  control,  C 

battle  management,  information  processing,  surveillance  sensors,  * 

!•>  intelligence  data  collection  and  handling,  solid  state  sciences,  s 

\  electromagnetics,  and  propagation,  and  electronic,  maintainability,  and  ? 

*  compatibility.  , 


